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ABSTRACT 

Progress in understanding the embedded stars in Lk Ha 225 has been hampered by their 
variability, making it hard to compare data taken at different times, and by the limited resolution 
of the available data, which cannot probe the small scales between the two stars. In an attempt 
to overcome these difficulties, we present new near-infrared data on this object taken using 
the ALFA adaptive optics system with the MPE 3D integral field spectrometer and the near- 
infrared camera Omega-Cass. The stars themselves have K-band spectra which are dominated 
by warm dust emission, analagous to class I— II for low mass YSOs, suggesting that the stars 
are in a phase where they are still accreting matter. On the other hand, the ridge of continuum 
emission between them is rather bluer, suggestive of extincted and/or scattered stellar light 
rather than direct dust emission. The compactness of the CO emission seen toward each star 
argues for accretion disks (which can also account for much of the K-band veiling) rather than 
a neutral wind. In contrast to other YSOs with CO emission, Lk Hoc 225 has no detectable Bry 
emission. Additionally there is no H2 detected on the northern star, although we do confirm 
that the strongest H2 emission is on the southern star, where we find it is excited primarily by 
thermal mechanisms. A second knot of H2 is observed to its northeast, with a velocity shift of 
-75 kms -1 and a higher fraction of non-thermal emission. This is discussed with reference to 
the H2O maser, the molecular outflow, and [Sll] emission observed between the stars. 



Subject headings: instrumentation: adaptive optics — binaries: general — stars: individual 
(Lk Hoc 225) — stars: variables: other — infrared: stars 
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1. Introduction 

Herbig Ae/Be stars (Herbig 1960) are commonly accepted as the intermediate mass (1.5-10 M ) pre- 
main sequence counterparts to T Tauri stars. Hillenbrand et al. (1992) separated Herbig Ae/Be stars into 
three distinct categories based on their spectral energy distributions, similar to T Tauri classifications. How- 
ever, the exact morphology that causes the infrared excess is under debate, with most models assuming an 
accretion disk within a large scale envelope (Natta et al. 1993; Hillenbrand et al. 1992; Di Francesco et al. 
1994). To further complicate matters, studies of these systems have shown that binaries appear to be as 
common in Herbig Ae/Be systems as in T Tauri systems (Leinert, Richichi, & Haas 1997; Pirzkal, Spillar, 
& Dyck 1997; Testi, Palla, & Natta 1998). Thus, coeval star formation of Herbig Ae/Be stars, or even Her- 
big Ae/Be stars with lower mass T Tauri stars, can provide a wealth of information about intermediate star 
formation and the role of these systems in clusters. 

BD+40 4124 (=V 1685 Cyg), one of the original Ae/Be stars of Herbig (1960), is at the center of 
a partially embedded cluster of at least 33 stars with infrared excesses (Strom et al. 1972). The region, 
part of the giant star forming region 2 Cyg at a distance of ~ lkpc, has a significantly younger stellar 
population than the nearby OB associations. One of the interesting aspects of this region is that both low 
and high mass stars are forming simultaneously, possibly induced from an external shock into the cloud 
core (Shevchenko, Ibragimov, & Chernysheva 1991; Hillenbrand et al. 1995). Although optical images of 
the field around BD +40° 4124 are dominated by the star itself, LkHa225 (=V 1318 Cyg) is arguably one 
of the more intriguing sources in the stellar cluster. The Lk Ha 225 system is a possible embedded binary 
with apparent separation of ~5000 AU oriented north-south. It has been proposed that the system consists 
of both a young low mass star, Lk Hoc 225 North, and an intermediate-mass very young Herbig Ae/Be star, 
Lk Hoc 225 South, based on luminosity arguments (Aspin et al. 1994). Many authors have suggested that 
the southern source is a center of star formation activity based on a variety of tracers. Aspin et al. (1994) 
presented a 800/jm map which shows a strong peak on this source, while Palla et al. (1995) found that both 
high density molecular gas traced by CSJ = 5-4 and the total gas column density traced by C 18 0/ = 2-1 
are centered in this region. More recently a 5" resolution CS / = 2-1 map by Safier & Looney (2001) shows 
that the densest gas is centred on the southern star of this pair. 

In this paper we present adaptive optics observations, of J, H, and K band imaging as well as integral 
field Kband spectroscopy toward Lk Hoc 225. With the highest spatial resolution to date in the near-infrared 
for this system we are able, for the first time, to make clear distinctions between the emission from the two 
stars. This allows us to highlight some intriguing differences in the K-band emission from the two stars, as 
well as a number of similarities. Moreover, we present the results of our unique ability to probe the inner 
binary region in a number of K-band emission lines. 

2. Observations & Data Reduction 

The observations described here were taken using the ALFA adaptive optics (AO) system (Kasper et al. 
2000) on the 3.5-m telescope at the German-Spanish Astronomical Centre on Calar Alto. While ALFA has 
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a laser guide star (Davies et al. 2000), here it was used in natural guide star mode. The wavefront reference 
was BD+40°4124 (niv=10.6, niK=5.6), which lies 36" from LkHoc225 and is the only star in the group 
brighter than my=12. While allowing a good correction on-axis, isoplanaticism at such a large angular 
distance is the main source of residual wavefront errors and severely limits the resolution achievable on our 
target. The only real alternative is to use Lk Hoc 224 which is much closer. However, while only 14" from 
the target it is too faint to correct on in the prevailing conditions; its magnitude has been seen to vary in 
the range 12.5 < V < 17.2 between 1985 and 1998 (Herbst & Shevchenko 1999). A discussion of the AO 
issues, with particular reference to the data obtained with Omega-Cass, is given in Davies et al. (1999). 

A K-band image of the field obtained with the near-infrared camera Omega-Cass in 1998, shown in 
Fig. 1, has a FWHM resolution of 0.25". This and the corresponding J- and H-band images (Fig. 2) were 
flux calibrated relative to star #3 in Hillenbrand et al. (1995). This star is the brightest unsaturated star in the 
image (K = 10.93 mag), but more importantly it is the least reddened with a J — K colour of only 0.22 mag, 
which tends to suggest it is the least susceptible to the variability associated with young stellar objects. 

We have also obtained integral field spectroscopy using 3D (Weitzel et al. 1996) during Oct 1999 at a 
pixel scale of 0.25", using AIM (Anders et al. 1998) as the interface between 3D and ALFA. The pixel scale, 
ambient seeing conditions (~1"), and distance to the reference star limited the final resolution to ^0.6". 3D 
is an integral field spectrograph that simultaneously obtains spectra for each of 256 spatial pixels covering a 
square field of view. The spectral wavelength range spans from 1.912/jm to 2.382,1/m at a spectral resolving 
power (R = X/AX) of 1000. Two pixel scales can be selected with AIM, 0"07 and 0"25. The latter was 
used for these observations, giving a field of 4" x 4". To observe both Lk Hoc 225 N and S, two pointings 
were necessary. The total on-source integration time for each was 600 s with individual frame integration 
times of 100 s. The same amount of time was spent off-source for sky background subtraction. The data were 
reduced using the set of 3D data analysis routines written for the GIPSY (Van der Hulst et al. 1992) package. 
This included wavelength calibration, spectral and spatial flat fielding, dead and hot pixel correction, and 
division by the reference stellar spectrum. Data cubes from the individual exposures were recentered and 
added using the centroid of the broadband emission. A mosaic of the two pointings was created using the 
relative coordinates of Lk Hoc 225 N and S given in Hillenbrand et al. (1995). Emission line and absorption 
line maps were extracted by subtracting a spline fit to the entire line-free continuum in each spatial pixel. 



3. Morphology & YSO class 

The system Lk Hoc 225 comprises the pair of stars in the south-east corner of Fig. 1, which are separated 
by 5", about 5000 AU at the distance of the cluster (980 pc, Hillenbrand et al. 1995). While this is rather 
large, nearly 5% of known YSOs do have such wide separations (Mathieu 1994), although it is not certain 
whether they are gravitationally bound or simply have formed close to each other from the same cloud core. 

The strong variability of this object, brightness changes of more than 2 mag in V and V-R colour 
changes of more than 2 mag within 1 year (Ibrahimov et al. 1988), means that it is not straight forward 
to compare observations at different times. Based on an R-band image, Aspin et al. (1994) identified 3 
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components to Lk Hoc 225 in a north-south line and a total separation of 5". The bright central component 
was also seen in the R-band image of Magakyan & Movsesyan (1997), distinct from both LkHoc225N 
and S. It is difficult to attribute this to a third star because the infrared appearance is quite different, with 
a long nebulous tail extending from the northern star. Fig. 2 shows that the tail is rather straight with a 
sharp kink at the north end, which may suggest it does not originate in the northern star but only appears 
to because it lies so close. The magnitudes of several points around Lk Hoc 225 (Table 1) and the colours 
marked in the JH-HK diagram indicate that not only is the whole ridge bluer than the two stars, but that the 
south end is bluer than the north end. If this is due to a decrease in extinction along the ridge of Ay=2 mag, it 
may explain why R-band images give the impression of a third discrete source. The colours along the ridge 
are easily achieved by reddening stellar light with a disk emission component by Ay~l 1. This is in contrast 
to the two stars which have a much redder H — K colour. Anomalous extinction laws have been observed 
towards a number of Herbig Ae/Be stars (Steenman & The 1989, 1991, and references therein) which can 
modify the direction of the reddening in the colour diagram, calculations by these authors have shown that 
this change is towards bluer H-K colours. Thus the position of the stars in the diagram appears to demand an 
additional significant (60-80%) dust contribution, which must arise in an envelope rather than a disk. Their 
position in the diagram suggests they are Class I young stellar objects (to which the Group I Herbig Ae/Be 
stars of Hillenbrand et al. (1992) appear to be analogous), or intermediate between Class I and II. Apart 
from Class sources (which do not emit at wavelengths shorter than 10/vm due to heavy extinction) these 
are the least evolved sources, with accretion disks, and surrounded by envelopes of gas and dust. The spectra 
presented in Fig. 3 tend to support this interpretation, having shapes suggestive of blackbody emission at 
~ 600 K, most similar to the Class I objects in the spectral atlas of Greene & Lada (1996). Such objects 
are associated with very high K-band veilings - defined as the ratio of excess to photospheric emission in 
the K-band, = F^ ex /F^. The excess emission is believed to come from hot dust or an accretion disk 
around the star and does not necessarily redden the star itself. Flares in such an accretion disk could explain 
the K-band variability of the two stars reported by various authors: the northern source was the brighter of 
the two in Sep 93 (Hillenbrand et al. 1995), Nov 93 (Aspin et al. 1994), and both Aug 98 and Oct 99 (this 
paper); while the southern star was brighter during Oct 91 and May 93 (Aspin et al. 1994). 

One possibility to explain the emission between the two stars is that the cloud core from which 
Lk Ha 225 formed condensed into a long filament. The continuum we see is tracing the remnants of this 
filament, arising either through direct emission from a string of young low-mass stars, or more likely since 
we do not see discrete sources, through light scattered from such stars. Polarisation studies would help to 
make progress on understanding this unusual morphology. 

4. H2 line emission 

The advantage of integral field spectroscopy is the ability to obtain line maps and spectra simulta- 
neously. Fig. 4 (centre) is a map of the H2 1-0 S(l) emission confirming what the spectra show, that the 
strongest emission is located directly on the southern star (we find no difference in their centroid positions), 
with a ratio S(l)/Bry>20. Additionally, there is an outflow/jet or extra source 1" to its north-east. For 
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reference, the left image (both greyscale and contours) is the integrated K-band continuum to indicate the 
positions of the 2 stars on the field. 

Aspin et al. (1994) observed weak emission from the northern star, while we have detected none. It is 
possible that a combination of seeing and large pixel size (3") meant that some of the flux from the southern 
star bled over to the pixels on which the northern star is located. With a finer pixel scale and high resolution 
our data do not suffer from this problem. 

A continuum-subtracted spectrum of the south star is presented in Fig. 5, in which 8 H2 lines are 
detectable at >3g, 4 from each of the v=l-0 and v=2-l series. The relative line fluxes are summarised 
in Table 2, and in Fig. 6 we have derived a relative population diagram for hot H2 molecules from the 
lines assuming the LTE value of 3 for the ortho/para ratio. While this assumption is not necessarily strictly 
correct, the diagram is a very useful presentation that gives an indication of whether the emission is thermal. 
We have further assumed that there is negligible differential extinction between 1.95-2.25/vm, in contrast 
to the A v =25±9mag estimated by Aspin et al. (1994) from the ratio of the lines 1-0 Q(3) to 1-0 S(l). 
We find that for Ay>9 implies the higher rotational levels have larger populations (even for the V = 1 
vibrational level which is derived from lines with high signal-to-noise). We suspect that Aspin et al. (1994) 
may have underestimated the continuum level longward of 2.35 fjm, as this is hard to determine at low 
spectral resolution due to the many H2 and CO emission lines: our spectra in Fig. 3 show that the continuum 
begins to rise rapidly at the long end of the K-band, the effect of which is to reduce the flux estimated in the 
Q-branch lines and thus also the derived extinction. 

From an analysis of H2 lines in an ISO spectrum, van der Ancker et al. (2000) concluded that there 
were two components: a C-shock (thermal component) with only mild (Ay~15) extinction, and a J-shock 
(non-thermal component) that is heavily (Ay~50) extincted. A reconciliation of these high extinctions with 
the low ones implied by our data may be possible if there is a significant extended H2 component, perhaps 
which we are beginning to see in Fig 4. 

Without any correction for extinction the v = 1 levels are well approximated by thermal excitation 
at a temperature of T ex = 1000-1100 K. Correcting for an assumed Ay=5mags of extinction raises the 
temperature derived from the v = 1 levels to nearly 2700 K. In order to make the most of even the weaker 
lines which have fluxes that are hard to measure we have taken the following approach. We have calculated 
the theoretical spectrum from 1.93-2.28/vm (to avoid the CO emission) using a combination of thermal 
(ortho/para ratio 3) and non-thermal components. For the latter, model 14 from Black & van Dishoeck 
(1987) (with ortho/para ratio ~1.7 and calculated for a density of 3 x 10 3 cm~ 3 and a incident UV radiation 
field of 3 x 10 3 times the background) was used. Since the relative line intensities are fairly insensitive to 
the density and strength of the radiation field, it provides a typical spectrum over the range of conditions 
that were considered; thus while we can estimate the fraction of H2 emission due to fluorescence, we cannot 
constrain the conditions for its excitation. Additionally, in our fitting procedure we have left the FWHM 
and velocity shift of the lines as free parameters. The results of the % 2 minimisation are given in Table 3. 
In summary, we find that the H2 lines at Lk Hoc 225 S are unresolved and have zero velocity shift; thermal 
excitation at 1070 K contributes 91% of the 1-0 S(l), and the remainder has a non-thermal origin. Two 
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temperature models are ruled out since a higher temperature component that is prevalent enough to account 
for the v = 2-1 lines would significantly modify the observed fluxes in some of the v = 1-0 lines. 

The calculation here assumes the low density regime, but on the other hand dense gas is a prerequi- 
site for H2O masers, one of which is seen in this region (Palla et al. 1995). The population diagram is 
indeed reminiscent of the high density photo-dissociation region (PDR) models of Sternberg & Dalgarno 
(1989), where H2 excitation occurs above the critical density at which collisional effects dominate the spec- 
trum. These authors found that for UV intensities more than 100 times the background and densities above 
10 5 cm~ 3 , the low vibrational levels become thermalised with a characteristic temperature of 1000 K. How- 
ever, their models showed that this extends to the v=2 level, exhibited as thermal ratios for the v = 2-1 as 
well as the v = 1-0 lines. This appears not to be the case for our data (Fig. 6), arguing against the scenario. 

We have extracted a spectrum of the H2 emitting region to the northeast of Lk Ha 225 S, in a 0.5" x 0.5" 
box. The same approach indicates that H2 lines in this region are resolved, with a quadrature corrected 
FWHM of 220km s" 1 , and are shifted by -75kms _1 . Only 81% of the 1-0 S(l) has a thermal origin, with 
a characteristic temperature of 980 K. That 19% of the 1-0 S(l) lines has a non-thermal origin implies that 
more than half of all the H2 cooling that can be deduced from these observations is through fluorescent 
emission. We cannot rule out the possibility that the broad line is due to two narrow components with 
different velocity shifts; but in any case there must be a significant component at -75kms _1 . This result is 
discussed in Section 6 with reference to the [Sll] line emission and the molecular outflow. 



5. CO emission 

The 2.3/mi CO 2-0 bandhead line map in Fig. 4 (right) clearly indicates that there is a compact spot of 
CO emission centered on each of the two stars and none elsewhere. In an analysis of the CO line emission 
seen with ISO, van der Ancker et al. (2000) used four rotational v = 0-0 transitions between 150-190,1/m, 
which have very low excitation temperatures and so trace cool gas heated thermally rather than fluorescently; 
hence they derived a thermal excitation temperature of only 300 K. If both the rotational and vibrational 
transitions of CO are excited collisionally (as argued by Carr 1989) then the CO 2-0 bandheads at 2.3/mi 
indicate an excitation temperature of more than 3000 K and a density greater than 10 10 cm~ 3 . 

Carr (1989), using their own and data from Geballe & Persson (1987), demonstrated that all the CO 
emitting sources then known also exhibited Bry emission, and that there was a correlation between them, the 
Bry luminosity being a factor of 2 or so greater than the CO 2-0 luminosity. Understanding this relationship 
has been difficult since the converse is not true, and only 25% of Bry emitters also exhibit CO emission. 
One unusual aspect of Lk Ha 225 S is that it breaks this correlation by exhibiting a ratio CO 2-0/Bry > 5 (we 
have not presented a Bry map since there is no detectable emission in 3D's entire field of view). Although 
Thompson (1985) proposed a model in which both CO and Bry emission originate in the same gas, Carr 
(1989) argue it is more likely that there are two components, hot (partially) ionised gas and cool (< 4000 K) 
neutral gas. Our observations can be explained if there is no hot component (or it is highly extincted), or 
if the temperature of Thompson's model is lowered. The concept of different lines arising in different gas 
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components can be extended to include the H2 emission in the southern star. The conditions needed to 
observe this line (1000 < T < 3000 K, density < 10 6 ~ 7 cm~ 3 ) are rather different to those for the CO, and 
so it must also arise in a different volume of gas. 

We cannot use the CO line profile from the observations presented here to distinguish between accretion 
disk or neutral wind models, but instead the CO morphology can help. The molecular outflow seen by Palla 
et al. (1995), for which they estimated a mass loss rate of 3 x 10~ 3 M yr _1 , would tend to support the wind 
model. To verify or veto this it is crucially important to know where the outflow originates. CO emission is 
seen only on the stars, inconsistent with an outflow that begins between them; indeed two distinct sources 
of CO emission is inconsistent with any model in which the CO arises in a single outflow. Additionally, the 
mass loss rate is also several orders of magnitude larger than models of winds predict from CO emission 
(Carr 1989). The alternative accretion disk model would naturally explain the compactness of the CO 
emission; and it also appears that it could simultaneously account for both the veiling and CO emission, 
at least in the southern source. In the northern source the CO flux is similar and could indicate similar 
conditions (density, temperature, accretion rate) in an accretion disk, but an additional warm dust envelope 
probably contributes to the higher veiling. 



6. Shocks, the Maser, & the Molecular Outflow 

From observations of the 12 C02-1 line at 1.3 mm, Palla et al. (1995) discovered a molecular outflow 
originating in Lk Hoc 225. They identified the origin of this with the H2O maser in this object which they 
confirmed to have a velocity shift of -80 km s -1 , and which appears to lie on the southern star. The coordi- 
nates given by Palla et al. (1995) and Hillenbrand et al. (1995) (B1950 and J2000 respectively) do not agree 
to better than a few arcsec, so we have used the 2MASS astrometry (accurate to <0.2") combined with our 
K-band image to determine accurate an position for LkHa225 S (a = 202030.59, 8 = 412126.0, J2000). 
The maser lies 0.4 ±0.2" north and 0.1 ±0.2" east of this star. If the transverse velocity of the maser is 
similar to its radial velocity, then it will not have moved by more than ~0.1" between the radio and our 
near-infrared observations. 

A number of different observations add to the complexity of the source. We have shown that there is a 
second weaker region of H2 emission to the northeast of Lk Hoc 225 S. It is also only 0.6" from the maser, 
and has a broader line width that is shifted by -75kms _1 , the same velocity as the maser. These are all 
suggestive of a link between the maser and H2 source. An important argument against the much brighter H2 
emission on the star itself being associated with the maser is that it shows no such velocity shift. A second 
line of reasoning is based on the maser model proposed by Elitzur, Hollenbach, & McKee (1989) in which 
H2O maser could occur in the ~400 K temperature plateau behind a J-shock. Elitzur (1995) further showed 
that C-shocks can also give rise to such masers, increasing the phase space to include preshock densities of 
10 6 -10 8 cm 3 and velocities 20-200 km s 1 . The spectrum of the H2 knot to the east of Lk Ha 225 S shows 
that more than half of the total H2 emission arises from molecules that have been non-thermally excited (a 
greater fraction than in the H2 emission coincident with the star). One way to achieve this is through H2 
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dissociation in a J-shock, if the shock velocity is more than 30^0 km s - although not faster than this 
since significant Bry emission must be avoided. It is the heat from re-formation of the molecules behind 
the shock front that produces the warm plateau ideal for maser formation (Elitzur, Hollenbach, & McKee 
1989). Based on these arguments, we suggest that this H2 knot and the H2O maser may be excited in the 
same shock system. 

The other emission which could be associated with the maser is the strong [Sll], which is also shifted 
by -80 km s -1 with respect to LkHoc225-N, and was observed between the two stars by Magakyan & Movs- 
esyan (1997). In this region they found almost no detectable Ha or [Nil] emission, putting tight constraints 
on the excitation temperature since the ionisation potential of hydrogen and nitrogen is 13.6 eV and 14.5 eV 
respectively, while for sulphur it is 10.4 eV. In their later observations (1995 rather than 1992) they also saw 
[Sll] had appeared at the location of the northern star. Such rapid variability in the emission lines effectively 
rules out shock-excitation since models show that the flow time through the hot shock structure is on the 
order of 100 yrs. We would argue that if this variability of the [Sll] is genuine then it is more likely to be 
in gas photo-ionised by one or both of the young stars (as discussed earlier, at least one may be late-B or 
early-A type). The -80 km s _1 velocity shift could be reproduced if the ionised gas was swept up in the 
outflow (while the average outflow velocity is 8kms _1 (Palla et al. 1995), it could be much higher nearer 
the stars). The variations might then be tied to changes of optical depth between the star and the gas, or be 
associated with flaring up of the accretion disk which would act as an intermittent source of higher energy 
photons necessary to produce [Sll]. 



7. Conclusions 

We have presented adaptive optics imaging and integral field spectroscopy of the deeply embedded and 
highly variable sources associated with Lk Ha 225 in the star forming complex around BD +40° 4124 . The 2 
bright stars have JHK colours and K-band spectra typical of Class I-II YSOs, with an additional component 
of emission coming from a more extended envelope rather than just a circumstellar disk. The extinction 
to the stars is uncertain, but seems to be Ay~ll mag. The ridge of continuum emission between the stars 
has bluer colours than the stars themselves, inconsistent with direct dust emission but appears instead to be 
stellar light that is extincted and perhaps also scattered. 

There are 2 sources of H2, one co-incident with the southern star, and one 1" to its north-east. For 
their line ratios to be physically meaningful neither can be reddened by more than Ay ~9 mags. In both 
cases most of the 1-0 S(l) emission has a thermal origin with a characteristic temperature of 1000 K, but the 
latter source does have a greater contribution from non-thermal emission. Its lines are broader and shifted 
by -75kms _1 , suggesting it is this knot which is associated most closely with the water maser (astrometric 
measurements place the maser almost equidistant from the star and this region). The variability reported 
by some authors in the [Sll] lines, which occur between the star and near the northern star, would seem to 
rule out shock excitation for this line. We suggest that it is a result of gas (which itself is caught up in the 
outflow) photo-ionised by one of the young stars, and that flares in the accretion disk or changes in optical 
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depth to the star cause the variations in line flux. 

Both stars show CO lines in emission, but in contrast to other CO emitting YSOs, there is no detectable 
Bry emission. The molecular outflow seen in the radio CO emission provides ample evidence for a neutral 
wind which could also give rise to 2.3/vm CO emission. However, two distinct 2.3/um CO sources is not 
consistent with a model of a single outflow; determining the exact origin of the outflow with respect to the 
K-band continuum emission would be a big step towards understanding the emission. On the other hand, an 
accretion disk model would naturally explain the compactness of the CO emission regions on each star, and 
it could account for at least some of the veiling seen. 

The authors are grateful to both the ALFA and 3D teams for their invaluable help, and the Calar Alto 
staff for their hosptality. We extend our thanks to Mansur Ibrahimov for providing a copy, and translation, of 
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Table 1 . Magnitudes of Lk Ha 225 



position 


magnitude 




K 


H 


J 


North star 


10.0 


12.7 


15.7 


South star 


11.5 


14.2 


16.6 


North end of ridge 


13.1 


14.9 


17.2 


South end of ridge 


13.5 


15.2 


17.2 



Table 2. Relative fluxes for observed H2 lines 



Line 


wavelength 


EJk* 


A u i a 


Relative Flux 




(pm) 


(K) 


(lO^s- 1 ) 


Southern star b 


Eastern emission c 


1-0 S(3) 


1.9576 


8365 


4.21 


5.30 


5.02 


1-0 S (2) 


2.0338 


7584 


3.98 


2.67 


3.53 


1-0S(1) 


2.1218 


6956 


3.47 


10.00 


10.00 


1-0 S(0) 


2.2235 


6471 


2.53 


1.93 


2.98 


2-1 S(l) 


2.2477 


12550 


4.98 


0.70 


1.35 


2-1 S(2) 


2.1542 


13150 


5.60 


0.26 


0.30 


2-1 S(3) 


2.0735 


13890 


5.77 


0.64 


0.96 


2-1 S(4) 


2.0041 


14764 


5.57 


0.38 


0.75 



a taken from the UKIRT Astro-Utilities page 
b ratios normalised to 1-0 S(l), with 1 a errors ~0.17 
c ratios normalised to 1-0 S(l), with 1 a errors ~0.12 
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Table 3. Derived parameters for H2 emission 



Quantity Southern star Eastern emission 



temperature of thermal emission 


1070 K 


980 K 


non-thermal fraction of 1-0S(1) 


9% 


19% 


non-thermal fraction of total H2 


38% 


59% 


FWHM of lines 


305 km s" 1 


375 km s- 1 


velocity shift 


-5kms _1 


^kms- 1 
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Fig. 1. — K-band image of the field around BD +40° 4124 , located top centre, plotted on a logarithmic scale 
(from Davies et al. 1999). The image has been smoothed with an adaptive filter so that bright sources are at 
full resolution while also allowing fainter sources to be seen. The source 5" northwest of BD+40° 4124 is 
a ghost not an astrophysical object. The other bright star is Lk Hoc 224, another Herbig Ae/Be star. The 
double system Lk Hoc 225 is the pair in the lower left of the image, with an obvious arm of continuum 
emission curling between them. These stars are 36" from BD +40° 4124 , which was used as the wavefront 
reference as it is by far the brightest source in the V- and R-bands. 
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Fig. 2. — Left - smoothed H-band image of the region around Lk Hoc 225 (H-band is used as it has the 
best contrast for the ridge). The ridge of continuum emission exhibits a sharp bend at its north end, which 
may suggest it does not in fact originate in the northern star. Right - colours of the 2 stars and 2 regions 
on the continuum ridge, as marked on the image, extracted in 0.8" apertures. The curved line starting near 
(0,0) shows the colours of dwarf stars, and the hashed region indicates the effect of adding models for disk 
emission; the reddening lines mark out the region where low and intermediate mass young stellar objects 
are expected to lie (adapted from Lada & Adams 1992). 
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Fig. 3. — K-band spectra of the 2 stars in LkHoc225, extracted in 0.88" boxes. The southern star exhibits 
strong H2 emission, while the northern star has none. The continuum in both cases is well matched by a 
450-500 K blackbody and a 20^10% contribution from an extincted A-type star. 
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Fig. 4. — Images of Lk Ha 225 observed with 3D and ALFA. Due to non-ideal observing conditions and the 
distance to the reference, the resolution achieved was ~0.6". Left: K-band (sum of all spectral channels) 
continuum showing faintly the emission curving to the west between the stars. Centre: H2 1-0 S(l) line map 
- there is strong emission both on the southern star and to its east, as well as some diffuse emission. Right: 
CO 2-0 bandhead 2.3/mi line emission map - the fluxes are similar, but the CO equivalent width for the 
northern star is much less. 
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Fig. 5. — Line spectrum of Lk Hoc 225 S after subtracting the continuum, with the features marked. 
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Fig. 6. — Population diagram for the hot H— 2 molecules in Lk Hoc 225 S, assuming a LTE ortho/para ratio 
of 3 and no extinction. While the assumption about ortho/para ratio is not necessarily correct, the diagram 
still shows clearly that the v = 1-0 transitions are dominated by thermal processes, while the v = 2-1 are 
predominantly non-thermal. The crosses are a best fit model to the spectrum, comprising a combination of 
thermal emission and UV fluorescence (assuming a density below 10 4 cm -3 ). The dashed line is the purely 
thermal part, which has a characteristic temperature of 1000 K. 



